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Abstract 

In this paper we have applied the cluster-expansion ansatz for the wave operator SI which in- 
corporates the orbital relaxation and correlation effects in an efficient manner. We have used both 
ordinary and normal ordered cluster operator (SI) to compute the ionization potential for different 
states of Ne, Ar, and Kr. The result in our calculation clearly shows how the relaxation and corre- 
lation effects play the role in determining the ionization potential, one of the interesting aspects in 
theoretical spectroscopy. 

PACS number(s). : 31.10.+Z, 31.15.Ar, 32.10.Hq, 31.15.Dv 

1 Introduction 

In spectroscopic studies ionization potential (IP), electron affinity (EA), excitation energy (EE) and 
double IP (DIP) are important in many aspects. For example, in all kind of electron spectroscopy, 
including photo-electron and Auger the IP, DIPs always play major roles. Specially, the knowledge 
of core IP provides insights into the choice and tunability of the ionizing beam in all kinds of electron 
scattering experiments like e — 2e, e — 3e, 7 — 2e, double Auger decay etc. These experiments are not only 
important to understand the physics but also to the development of new fields like quantum computing 
thorough the idea of entanglement involved in the process. Hence accurate theoretical determination 
of IP for core as well as valance electrons may give some impetus to the physics underlying in those 
experiments. 

Open-shell Coupled Cluster (OSCC) method is used since several years to compute the IP, EA p]|2lEl 
and DIPs PIE] - usually employing complete model space (CMS) [13 IE] In 1980's coupled cluster (CC) 
methods have been generalized to handle the incomplete model space (IMS) [21E1C1IH] m a size-extensive 
manner. In this paper we have tried to demonstrate the results arising from the algebraic structure of 
the expansion of the wave operator Q = exp(T + S). We have used both the ordinary and normal ordered 
cluster expansion of the wave operator f2. Normally the operator T corresponds to core-correlation. 

Lindgren jH] has shown how the normal ordered cluster expansion of the wave operator Q can include 
the correlation effects. Here we will show in a sophisticated way how the relaxation effects, the many- 
body effect of different kind can also be taken into account in the calculation of ionization potential 
within the algebraic structure of the semi-normal ordered cluster expansion of the wave operator fl. The 
form of the cluster expansion takes care the relaxation effect in a compact and efficient manner and it can 
be easily demonstrated how the coupling of different many body effects like relaxation and correlation 
effects come into the picture automatically. In an earlier publication Chaudhuri et. al. [X(]] have explored 
this techniques and applied to molecules. 

In this paper we have calculated ionization potentials (IPs) of different states of Ne, Ar and Kr. IPs 
are calculated for core and valance orbitals of those closed shell atoms using two different approaches. 
Section [2] deals with the theoretical formulation of the two different methods of calculation using the 
exp(S) and {exp(S')} theories respectively, where {• • •} denotes normal ordering. Section presents 
the results obtained from our calculation and a detailed discussion about the physical aspects of the 
two theories. In section 0] we have made the concluding remarks and the important findings of our 
calculations. 
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2 Theory 



It is a well established fact that the orbital relaxation effects are more important than the correlation 
effects in the core hole ionization (CHI) process. Since in CHI the relaxation effects are stronger than 
the correlation effects, perturbative treatment will not be a efficient one because in higher orders it will 
be very cumbersome and that is the reason to choose the non-perturbative approach like CC. 

Let us start from the TV-electron Dirac Fock (DF) determinant | <&o ) - The (N— 1) ionized determinants 
|$ Q ) are obtained from the DF reference state where one electron is ionized from the orbital a. Here the 
labels a, /3 . . . denote holes and p,q . . . denote particles. According to Thouless[ll], the orbital relaxation 
corrections can be introduced through an exponential transformation involving one body operator of the 
form: 

|$r d ) = exp(S?)|$ a ) = exp(S?)a a |*o>, (1) 

where 

S? = E (pN Q |7)^a 7 , (2) 

76|*a) 

provided 

($r d i$«> ¥= o. (3) 

The amplitudes (p\s a \j) can be uniquely determined and orbital relaxation effects may be induced 
through the above exponential transformation. Like Eq.JQl, \H>), the combination of the determinantal 
states |$ a ) also transforms in a similar way and we have the new function having modified determinant 
with the associated modified orbital set for each determinant. We write it as 

|$r d ) = £cxp(S?)P Q |*> =^C a expK)|$ a ) I (4) 

a a 

with 

|*}=EC7 Q |$ Q ), (5) 

a 

where C Q 's are the coefficients of \& and the P a 's are the projectors onto the <S> a . 

Now our aim is to define an expansion which is analogous to Eq.Q. To do so, we introduce a 
valence universal wave operator Q, |9; for both neutral and the (N-l) electron ionized model space states, 
satisfying the Fock-space Bloch equations El Ej 

HQP (k) = nP {k) H ef fP (k) . (6) 

where and PW are the projectors onto $o and {$q}, respectively. We invoke a normal ordered 
exponential ansatz [2j for Q as 

Q = Q c n v = exp(T) {exp(S*)} (7) 

that satisfies 

|$)=^fiC Q |$ a ). (8) 

a 

The cluster operators T and S appearing in Eq.Q are defined as 

r = X>l*ih r >+ E (pqlhhP) {alalapa 7 } , (9) 

and 

s = £>|*?|7) + E <H4 It/3) {44 W ■ ( 10 ) 

Here {■ • ■} denotes normal ordering with respect to <f>o. Due to the normal ordering, S acting on <£> Q will 
give zero unless the determinantal state has a hole in a. For convenience, we express the cluster operator 
S in terms of correlation (S c ) and relaxation (S r ) operator as 

S = S r + S c (11) 

where 

S r = ]TH s ?l7> + E <H fi 2 17«) {44a«a-y} , (12) 
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and 

S c = ]T {pal^lT/SXototo^}. (13) 

Figures ^)-lQb) represent the Hugenholtz diagrams corresponding to the terms given in En. (T2)l and 
figure (2b) represents En. (T3j) for (3 ^ a. 

Since {S} 2 and higher powers of S contains more than one valence destruction operator and 
contraction among S operators is not allowed, Eq.Q reduces to 

n = n c n v = n c {i + s r + s c }. (14) 

It is evident from En. (141 that {exp(S')} = {1 + S} can not induce relaxation corrections via Thouless- 
type exponentiation. 

Let us now consider expansion of fl where f2 is an ordinary exponential ansatz, i.e., 

Q = fl c n v = exp(T) exp(S) (15) 
Substituting Eq.Jll} in Eq.JEJl, we get 

n = n c [i + {s c } + {s r } + j {IFW} + i {s r s r s r } + j {WWs } + ■■•]. (16) 

This expansion looks similar to the normal ordered expansion with the extra property that the con- 
tractions like S r S r involving the orbital a is allowed. It can easily shown that in the above expression 
(EqEUi terms like {S r } 2 , {S c } 2 etc. do not contribute because these terms contain more than one 
valence destruction operator a* . Introducing an operator a r through the relation 



a 

1 n 

n—l 



J2^{S r ---S r ---S r } , (17) 



we can rewrite Ea. ljlfifl in normal order as 

n = n c {exp(s c + a r )} = n c n v . (is) 

En.(T8| clearly indicates that the ordinary exponential ansatz exp(5) is capable of inducing the or- 
bital relaxation effects through Thouless-type exponentiation. The relaxation effects arising from the 
contraction of terms like S r S r are shown in figure [2] where in (2K) (3,5 ^ a and in figure (2b) P, 7, 5 ^ a. 

To compute the energy difference, we insert Eo. (14t /Ea. tlSjl into Eq.© and proceed hierarchically 
from k = to k = 1. For k = 0, Eq.® reduces to closed shell CC equation for core-cluster amplitudes 
T. For k = 1, we premultiply Eq.® by O^T 1 (= exp(— T)), introduce a connected dressed operator 
H = fl~ 1 Hil c — E gr and obtain the Bloch equation for energy difference 

Hn v pW = n v p {1) H eff P il \ (19) 

where 

Hef/ = H eff - E gr . (20) 



3 Results and discussions 

The results obtained by our calculation of IPs of Ne, Ar and Kr using two different approaches are 
presented in table Figures 03 El and show the relative error for all the states of Ne, Ar and Kr 
respectively, where as the figure shows the relative error in IP for Is states explicitly of the systems 
studied. From tableland figures Q] and El it is clear that the agreement of IPs obtained by the exp(5) 
theory is better than {exp(S')} theory for the core states. This nature can be interpreted in terms of 
the different many-body effects, namely relaxation and correlation effects. For deep lying cores, the 
relaxation effects are stronger than the correlation effects, which is taken care of by the exp(S) theory as 
compared to the {exp(S)} theory (FignjGJ. Where as for high-lying core orbitals the relaxation effects 
are comparatively less than the deep core orbitals like lsi/ 2 and 2si/ 2 - And thus the results obtained by 
the {exp(S')} theory performs better than the exp(S) theory for high lying states like 3p and ip. These 
can be interpreted in a better way if we look at the Thouless's treatment JH- When the relaxation is very 
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(a) (b) 
Figure 2: Diagrammatic representation of the cluster operator S 



Table 1: Ionization potentials (in eV) of Ne, Ar and Kr atoms. Entries within parentheses are from 
UCC(3). 
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strong, a perturbative treatment is very cumbersome in higher orders and also inefficient. The widely 
used normal ordered cluster expansion approach is not compact enough to include the relaxation effects 
in an adequate manner in the theories to calculate IP, EA, EE and DIP. To overcome this problem we 
have explored the generalized Thouless's treatment and have calculated the results by using the exp(5) 
theory. 

In table^the results in the parentheses are obtained from UCC(3) and it is clear that for these kind 
of problem the introduction of the unitary excitation operator does not improve the results much. The 
details and the results using UCC(3) |TIJ will be communicated shortly. 

If we minutely observe the results for Kr we can see that the relaxation and correlation effects are 
canceled out for almost all the states. In other words, these two effects are compensated by each other. 
This becomes more clear if we see the Koopman's energy for this system given in tabled 

Besides the experimental values, IPs from <5-SCF are also listed in table The 5-SCF IPs are 
computed using GRASP programme In this procedure, the electron binding energy corresponds to 
the ath orbital is obtained from 

IP{a)=E(N)-E a (N-l), (21) 

where E(N) and E a {N — 1) are the SCF energies of N and (N — l)-electron systems respectively. The 
present calculation shows that 5-SCF estimates the core IP more accurately than the valence IP. This 
is in accordance with our expectation. Note that the differential electron correlation which is important 
for valence IP can not be incorporated through (5-SCF procedure. Since the orbital relaxation effect can 
be introduced through 5-SCF procedure, this yield reasonably accurate core IPs. 



4 Concluding remarks 

In this paper we have presented the calculation of ionization potentials for core and valance orbitals of 
some rare gas atoms, namely Ne, Ar and Kr. We have explored both the exp(S') and {exp(5')} theory 
to obtain the results. Along with the normal cluster expansion ansatz we have also obtained the results 
by using unitary cluster operator theory (UCC) for both the formalisms. 
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5 SCF exp(S) UCC-exp(S) {exp(S)} UCC-{exp(S)} 

Figure 3: % error for different states in computing IP of Ne using different methods 




8 SCF exp(S) UCC-exp(S) {exp(S)} UCC-{exp(S)} 

Figure 4: % error for different states in computing IP of Ar using different methods 
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Figure 5: % error for different states in computing IP of Kr using different methods 
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Figure 6: % error for different Is states in computing IP of Ne, Ar and Kr using different methods 
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The close comparison between the results obtained by the two different approaches are explored 
and discussed in details. The figures containing the computational errors for different theories help us 
to understand the trends of the obtained results. This calculation clearly demonstrate the importance 
of different many-body effects, namely differential electron correlation and orbital relaxation effects in 
determining ionization potential of closed shell atoms. 
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